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. ABSTRACT 

o: 

, Using a sample of 123 X-ray clusters and groups drawn from the XMM-Cluster Survey 

first data release, we investigate the interplay between the brightest cluster galaxy (BCG), its 
black hole, and the intra-cluster/group medium (ICM). It appears that for groups and clusters 
with a BCG likely to host significant AGN feedback, gas cooling dominates in those with 

■ Tx > 2 keV while AGN feedback dominates below. This may be understood through the 
\—i [ sub-unity exponent found in the scaling relation we derive between the BCG mass and cluster 
5^ . mass over the halo mass range lO^'' < M500 < lO^^M© and the lack of correlation between 

radio luminosity and cluster mass, such that BCG AGN in groups can have relatively more 
energetic influence on the ICM. The Lx — Tx relation for systems with the most massive 
BCGs, or those with BCGs co-located with the peak of the ICM emission, is steeper than that 
for those with the least massive and most offset, which instead follows self-similarity. This 
is evidence that a combination of central gas cooling and powerful, well fuelled AGN causes 
the departure of the ICM from pure gravitational heating, with the steepened relation crossing 
self-similarity at Tx = 2 keV. Importantly, regardless of their black hole mass, BCGs are 
more likely to host radio-loud AGN if they are in a massive cluster (Tx > 2 kcV) and again 
co-located with an effective fuel supply of dense, cooling gas. This demonstrates that the most 
massive black holes appear to know more about their host cluster than they do about their host 
galaxy. The results lead us to propose a physically motivated, empirical definition of 'cluster' 
and 'group', delineated at 2keV. 

Key words: galaxies: clusters: intracluster medium - galaxies: active - galaxies: elliptical 
and lenticular, cD 



1 INTRODUCTION 

* E-mail: j.p.stott@durham.ac.uk Galaxy clusters and groups are important probes of cosmology as 



2 J. P. Stott et al. 



scendants of the first over-densities to collapse after the big bang. 
They can be parameterised by the X-ray emission properties of their 
hot intra-cluster medium (ICM), with both the luminosity of the 
X-ray emission (Lx) and the temperat ure of the gas (T x) scaling 
with the total mass of the system (e.g. IPratt et al.ll2009[) . Lx also 
scales with the integral of the square of the gas density while Tx 
indicates the virial temperature of the system. For a self-similar 
relation, wh ere grav itational heating is the only source of energy, 
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2-3 (lMushotzkYlll984|:lEdge & Ste wart 1 99 ll:lMarkevitctill 19981 : 
lAmaud & Evrardll 19991 : Iprattet al.ll2009i : iMittal et alj|20 lib which 



demonstrates that s econdary and non-gravitational effects are in- 
fluencing the ICM ( Allen & Fabian 1998: Muanwong et al. '200^; 
Randall et aTllIoOZ: .Rowlev et aL.2004. : .Magliocchetti & BruggenI 
2007l:lHartlevetalj2008l) . 



As both Lx and Tx are used to calculate cluster mass but 
probe different physical properties of the gas, it is possible to as- 
sess the state of the ICM by studying the relation between them for 
a sample of clusters. It is important to understand this Lx — Tx re- 
lation and any definable sub-samples within it, as X-ray mass scal- 
ing relations are used to derive cluster masses and are employed 
as complementary probes o f cosm ology to th ose derived from 
other methods ( Borgani e t al.ll200ll:[s tanek et al]|2006l : lAllen et al.l 
I2OO8I : IVikhlinin et al.„2009l: ISahlen et al..,2009ll . A further motiva- 
tion is to explore the physical processes that act within the ICM 
itself. For example, we seek to understand the process that pro- 
vides an entropy excess to stop the gas in the centres of clus- 
ters from ca tastrophically over-cooling via the classical cooling 
flow model ( lFabian|[T994h . The leading candida te for this is the 
energy output by active galactic nuclei (AG N, Tucker & David 
1 19971: iMcNamara & NulsenI I2OO7I: lOitti et aljliol l|). This is also 
often invoked in one form or another to reconcile the dark mat- 
ter mass function with the galaxy luminosit y function jBest et al] 
l2006l : ICroton et al.ll200d : lBower et al.l|200 d) and is thus a cosmo- 
logically important energy source across a range of mass scales. 

In terms of gas cooling, for dynamically relaxed clusters that 
have a cool-core present (i.e. the cooling time is significantly 
less than a Hubble tim e) the central density of gas and thus Lx 
will be relatively high (Fabian '1994'). This manifests itself in part 
through an offset or steepening of the Lx — T x relation for cool- 
core clusters rela tive to the whole population dFabian et al .11 19941 : 
iMittal et al.l201 ll) . However when corrected using both X-ray spec- 
tra and imaging or ignoring the central re gions, the self-similar 
result can be r ecovered for some clusters 1 Allen & Fabian|[l998l : 
iMaughan et al.ll201 1). In contrast, unvirialised systems, for exam- 
ple those undergoing a cluster-cluster merger, will have disturbed 
X-ray emis sion and a Tx that is less indicative of th e total mass of 
the system jRandall et al.l2002l : lRowlev et al.ll2004l) . 

As mentioned above, the galaxies themselves can influence 
the gas through non-gravitational processes due to the presence of 
AGN which can heat and even blow the gas out of the central re- 
gions of clusters and groups, and their progenitors. The effect of 
radio-loud AGN activity on the ICM has been seen in a number of 
detailed X-ray and radio observations of local clusters and groups 
where the radio lobes of the AGN are found to inflate large cavi- 
ties i n the X-ray emitting gas, demonstrating t he power of this pro- 
cess dChurazov et al.ll200ll : lBirzan et al.ll2004b . In the cluster mass 
regime the power in these AGN cavities matches the luminosity of 
the central cooling ICM, whereas in low mass clusters and groups 
the cavity power is found to exceed that of the ICM, a nd thus the 
energy injection from AGN can dominate in the core dCitti et al] 



sity of clusters, and therefo re Lx, while the Tx is less affected, 
as is t he case in simulations dPuchwein et al.l2008l : lMcCarthv et al] 
'2OIOI). Alternatively this can be thought of as heatin g the gas for 
a giv en Lx as has been reported observationally (Croston et al.l 
I2OO5I) . 

There is further observational evidence connecting the galaxy 
scale to the cluster scale, as the broad-band optical and near- 
infrared properties of the central galaxies in groups and clusters 
(hereafter BCGs) have been found to correlate with the X-ray lu- 
minosity (Lx ) and X-ray temperatur e (Tx ) of the hot ICM of their 
hosts, albeit with significant scatter ('Edge"l99l': 'Collins & MamJ 
1998; Lin & Mohr 2004: Popesso et al. 2007; Whiley et al. 200l; 
Stott et al.ll2008l : iBrough et alj|2008l : IMittal et al.ll2009 ). As these 



X-ray properties correlate to first order with the total mass of the 
system, BCG stellar mass must correlate with the dark matter halo 
mass of the cluster. This may be indicative of the hierarchical 
build up of the galaxies and their host halos over cosmic time, 
as these masses are also f ound to correlate in theoretical models 
dPe Lucia & Blaizotll2007l) . 

We may therefore ask, are the BCGs influenced by and 
influencing the ICM, and if so how? We can infer from the 
AfeH — cr relation, in which black hole mas s is found to corre- 
late with the host galaxy ve locity dispersion dFerrarese & Merrit^ 
|200(]| ; iGebhardt et alj|2000l) . and the relationshi p between black 
hole mass and galaxy bulge mass (A/rh — AfR,,i^n. lMagorrian et al.l 
Il998l) . that all BCGs contain supermassive black holes. How- 
ever because of the duty-cycle of supermassive black hole ac- 
cretion, only a certain fraction are AGN at any one instant. Ob- 
servationally, cluster galaxies and BCGs in particular are more 
likely to host AGN and these are more likely to be powerful ra- 
dio sources, having the highest duty-cycles of any class of galaxy , 
with duty-cycle correlating with stellar mass dLin & Mohj|2007l: 
iBest et al.ll2007l) . Clustering and environment studies of radio-loud 
AGN consistently find that they typically reside in halos of mass 
> IO ^^Mq (e.g.. Best 2004; Wake et al. 2008; Mandelbau m et al.l 
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1201 ID . AGN processes such as these can lower the central gas den- 



l2009l ; iHickox et alj|2009l : IPonoso et al.ll2010l) . precisely the halos 
with virialised hot atmospheres for which AGN heating is required 
in cosmological models. There is also a correlation between BCG 
stellar mass and the maximum radio lu minosity observed for that 
mass of galaxy, albeit with large scatter dLin & Mohr 2007.) . This 
can be understood in terms of the AfsH — Meuige relation, and 
the fact that AGN power output correlates with black hole accre- 
tion rate, which i s a fun ction of black hole mass (e.g. Eddington or 
Bondi accretion, Bondi 1952). The significant scatter is due to ob- 
serving the AGN at different stages of the accretion cycle between 
observationally quiescent and radio-loud AGN. One might there- 
fore expect that the most massive BCGs will have supermassive 
black holes able to inject the most energy into their surroundings 
and thus most strongly influence the ICM. 

This may imply that clusters and groups containing the most 
massive central galaxies will have a lower Lx for a given Tx , as 
they will be those that possess black holes able to inject the most 
energy into the ICM. However, the process of AGN accretion and 
then energy injection into its surroundings is often described as 
a feedback process that regulates its own accretion rate by heat- 
ing and/or mechanically removing i ts own fuelling gas supply 
dBest et al.' 2006; Bo wer et alj|2006l : iMcNamara & NulserJ I2OO7I: 
[Booth & Schaye 201(J). In-falling gas fuels the accreting black 
hole, producing radiation driven winds from the accretion disk or 
mechanical jets of radio emitting material, that act to stop the gas 
from the surrounding medium from cooling onto the accretion disk. 
This requires the energy output to be coupled in some way to the 
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gas cooling rate. So well-fuelled AGN should be more powerful 
than those with less available gas, but these well-fuelled AGN will 
necessarily be in higher gas density environments, potentially set- 
ting up a balance between AGN feedback and gas cooling s ufficient 
to sto p the over-cooling predicted from radiative processes jpabianl 
Evidence for the co-location of a central AGN with a cold 



gas fuel supply comes from CO detections in BCGs at the cent res of 
cool- cores that display both radio and optical line emission teded 
l200lh . 

Previous work has provided tantalising clues on the inter- 
relation between the ICM and the galaxy population. For exam- 
ple looking at the Lx — Tx relation for clusters that contain 
radio emission, those that have extended radio structure have a 
significantly steeper relationshi p than those that are unresolved 
jMagliocchetti & Bruggenll2007 ). In t erms of cool-core clusters, as 
mentioned above, Mittal et al. 201 ih find that those with strong 
cool-cores have a significantly steeper Lx — Tx than those lacking 
a cool-core and that the presence of a cool-core indicates a radio- 
loud central galaxy. There is clearly a relationship between both of 
these observations, in that it seems that to guarantee a central AGN 
requires a strong cool-core, and cool-cores and/or AGN lead to a 
steepening of the Lx~Tx relation. This demonstrates perhaps that 
cool-cores can dominate over AGN feedback in the cluster envi- 
ronment, while AGN become more significant in the group regime . 
Further evidence for this picture is presented in lLin & MohJ bOOTi) 
who find that the ratio of energy injected by AGN to the thermal en- 
ergy of the ICM decreases with increasing cluster mass and recent 
observations of the radio heating of the X-r ay gas around individual 
early-type galaxies jPanielson et al.ll2012h . 

In this paper we perform the first comprehensive study of 
the interplay between the stellar population of the BCG, its super- 
massive black hole and the ICM by investigating the optical and 
radio properties of BCGs within the context of the cluster/group 
Lx — Tx relation. The sample used consists of 1 23 z < 0.3 
clusters taken from the XMM C luste r Survey (XCS, iRomer et al] 
l200lHLlovd-Davies et al.l201lHMehrtens et ailioilh . a serendip- 
itous survey of clusters from the XMM-Newton archive, a large 
sample which contains a range of systems from the group scale up 
to massive clusters, which is key for this study. The combination of 
both optical and radio observations means we can assess both the 
time averaged (through the assumption of the A/bh — A/euigc re- 
lation) and instantaneous effect of AGN feedback on the ICM. The 
results of this paper are discussed assuming that all BCGs contain 
supermassive black holes, obeying the A/bh — A/euigc relation, 
that are potential AGN but whose duty-cycle means that we can 
observe only a fraction of them to be accreting at any one time. 
Where appropriate we compare our results to those from the Over- 
whelmingly L arge Simulations (OWLS) run which inc orporates 
AGN feedback jSchave et alJl2010l : lMcCarthv et al.ll201(]|) . 

We begin the analysis in ^3. II with BCG - cluster X-ray scal- 
ing relations, in order to derive a BCG mass to host dark matter 
halo mass relation. With the assumption of the A^bh — A4"Buigc re- 
lation, this allows us to assess the relative significance of the BCG 
black hole to that of the entire system. In i]3.2| we study the effect 
of BCG mass, and thus the time averaged effect of its black hole, 
on the ICM through the location of the cluster within the Lx — Tx 
relation. Under the assumption that black holes also require an ef- 
fective supply of fuel for accretion, in i]3.3| we also investigate the 
effect of the projected distance of the BCG from the peak surface 
brightness of the ICM on the Lx — Tx relation. We then move on 
to investigate the instantaneous effect of AGN in i]3.4| to look for 
the key factors affecting the BCG radio-loud fraction and the effect 



of radio-loudness on the Lx — Tx relation. Finally, in §4T|we use 
simple energetic arguments to explain our findings. 

A Lambda Cold Dark Matter (ACDM) cosmology (fim = 
0.27, = 0.73, Ho = 70 kms"^ Mpc"^) is used throughout 
this work. 



2 SAMPLE AND DATA 

The galaxy cluster sample and the ass ociated X-ray propert ies are 
taken from the XCS first data release (Mehrtens et al. '201 1*), with 
the data analysis techniques described in Llovd-Davies et al. ( 201 1]) 
and we therefore point the reader to those papers for the sam- 
ple details. In order to compare the optical and radio properties 
of these systems, the additional selection criteria are that the clus- 
ters lie within the Sloan Digital Sky Survey Data Release 8 (SDSS 
DR8) footprint and, where appropriate in this paper, also within 
the Faint Imag es of the Radio Sk y at Twenty centimetres survey 
(FIRST survev. lBecker et aljl995h . which has a footprint very sim- 
ilar to SDSS. A redshift cut of z < 0.3 is also imposed to ensure 
good optical photometry and completeness of radio detections. The 
former XCS-SDSS sample consists of 123 clusters and the latter 
XCS-SDSS-FIRST sub-sample consists of 103 clusters. The XCS 
has a complicated selection function because of its serendipitous 
nature and heterogeneous data depths, but ¥igme\I\(upper panel) 
demonstrates there is no evidence for a strong trend between Lx 
and redshift due to a flux limit and, importantly, that we find a sim- 
ilar dynamic range in Lx over the entire redshift range considered 
here. 

The key advantages of the XCS first data release are the 
sample size and its sensitivity, providing a large range in Lx 
(0.007 — 14.0 X 10**ergs~^, for our sub-sample) and thus an 
excellent dynamic range in a property related to gas density, es- 
sential to probe the effects of AGN on the ICM. The wide range 
in Tx (0.3 — lOkeV, for our sub-sample), and therefore mass, al- 
lows us to analyse a significant number of group scale systems, 
where AGN become more energetically important. Importantly, 
for a study of clusters that may also harbour X-ray bright AGN, 
the X-ray data have had any detectable point sources removed, as 
they were detected within, or in the vicinity of, ~ 25% of the 
z < 0.3 sample. These are dealt wi t h usin g the XAPA vision 
model described in iLlovd -Davies et al. l|), that can detect and 
remove point sources embedded in extended sources and differen- 
tiate between blended point sources and genuine extended sources. 
The X-ray properties and their associated errors can be found in 
iMehrtens et alj ( I2OI ih but we illustrate the X-ray errors in the Ap- 
pendix Figure lATI The X-ray luminosities (Lx) quoted throughout 
the paper are bolometric and are calculated within -R500, the ra- 
dius at which the cluster mean gravitational mass density drops to 
500 X that of the critical density. The temperatures are also calcu- 
lated within R500, using spectral fitting c arried out w ith XSPEC 
and the maximum likelihood Cash statistic l lCashll 19791) . The mod- 
els used i n the fitting include a photoelectric absorption component 
(WABS: [Morrison & McCammonI Il983h to sii nulate the nH ab- 
sorption and a hot plasma component (MEKAL: lMewe et aljl98^ 
to simulate the X-ray em ission from the ICM (for full details see 
iLlovd-Davies et al.ll201 ih . 

We calculate M500 masses for the clusters using an evolv- 
ing mass - Tx scaling relation as described in Sahlen et al. ( 200' 
and the virial mass-to-Afsoo conversion from iHu & Kravtso 
(2003), with the average conversion from virial mass to Afsoo 
being a factor of 0.79 ± 0.02. The M - Tx relation is nor- 
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malised using HIFLUGCS data for the local cluster population 
teeiprich & B6hringej|2002h . Implicit in the mass determination 
is also the assumption of an NFW halo profile with a concentration 
parameter of 5 ( Navarr o et al. 1995). As the redshift range is rela- 
tively small, the result of this derived M500 — Tx scaling is well fit 
by the relation: 

log (M500) = 1.432(±0.011) log {Tx) + 13.722(±0.005), (1) 

with M500 in Mq and Tx in keV. The XCS sub-sample in this 
paper covers the mass range 10^^ < M500 < 10^^ . 

The BCG optical properties are extracted from the SDSS 
DR8 i-band model magnitude photometry which assumes a de 
Vaucouleurs profile fit to obtain a total magnitude. Such a fit is 
appropriate for the BCGs in our sample as we find the median 
Sersic n = 4.0 ± 0.4 for a sample of X-ray selected clusters at 
z ~ 0.2 observed with the Hubble Space Telescope jStott et aEI 
I2OI ih . The magnitudes are converted to a lumin osity (Lj-spss) us- 
ing k and passive evolution corrections from a lBruzual & Charloj 
( I2OO3I) simple stellar population (SSP) with a Chabrier IMF, a for- 
mation red shift, Zf = 3 a nd solar metallicity, a model appropriate 
for BCGs ( Stott etal. 2010). A correction for Galactic extinction 
using the dust maps of Schlegel et al. ( 1998) is used. We also calcu- 
late a stellar mass for the BCGs using an i-band mass-to-light ratio 
of 1.91 derived from the same SSP. Figure[T](mWii/e panel) demon- 
strates that optical incompleteness is not an issue out to z = 0.3. 

The radio data are taken from the FIRST survey catalogue, 
which has a detection threshold of 1 mJy. From catalogue match- 
ing and visual inspection of the FIRST imaging we obtain the flux 
from the radio source centred on the BCG and also sum the flux 
from all radio sources associated with the BCG (i.e. central point 
source and extended lobes). We note that all of our radio sources 
have more than just a central component and can therefore all be 
classified as extended. The 1 .4 GHz radio luminosity of the central 
and total components (I/Radio-ccn and I/Radio-tot) of the BCGs 
is calculated with a k correction assuming a spect ral slope where 
flux, oc and a = 0.72 l lCoble et al.ll2007h . To ensure the 
radio detected fraction of the XCS-SDSS-FIRST sample is com- 
plete to z — 0.3, we plot I/Radio-ccn against redshift in Fig- 
ure [T] (lower panel) and implement a luminosity limit such that 
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nature of the sample. Hereafter I/Radio-ccn > 2 x lO^'^WHz^^ 
defines the 'radio-loud' sample, with the remaining XCS-SDSS- 
FIRST clusters being 'r adio-quiet' . We n ote, t hat this is similar to 
the lower limit used in Best et al. ( 2007) and iLin & Mohil j2007l) 
(1 X 10'^'^ W Hz^^) with the former finding that the contamination 
of the radio-loud BCG sample with star forming galaxies is min- 
imised to < 1% using this criterion. The most luminous source 
in the sample is the cluster associated with the well-studied ra- 
dio galaxy 4C -1-55.16 with LRadio-tot = 1.4 x lO^^WHz"^ at 
z = 0.24. 

The BCG selection for a cluster is usually obvious from visual 
inspection of the images as they are the prominent galaxy closest 
to the X-ray centroid and often have a cD-like profile. However we 
choose to formalise this by studying the tip of the red sequence in 
the colour magnitude relation. For each cluster we identified the 
red sequence with g — r colour and selected the brightest galaxy 
from the r-band magnitudes of all the red sequence galaxies within 
a projected distance of 500 kpc from the cluster X-ray centroid as 
for approxim ately 95% of clusters the BCG lies within this radius 
dLin & Mohj|2b04) . 

A table of optical and radio data from this paper can be found 
at: \http://www.astro.dur.ac.uk/^jps/Stott2012/Stott2012.cat\ 
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Figure 1. Upper: The cluster/group bolometric X-ray luminosity plotted 
against redshift for the XCS-SDSS sample. The radio-loud BCGs ai'e plot- 
ted as open diamonds. This demonstrates that Lx does not strongly depend 
on redshift. Middle: BCG SDSS i-band luminosity plotted against redshift. 
The radio-loud BCGs are plotted as open diamonds. The sample shows 
no X/i— SDSS - redshift dependence. Lower: BCG central radio luminosity 
(I/Radio-ccn) plotted against redshift for all radio detections in the sample. 
The horizontal line at iRadio-ccn = 2 X lO^^W Hz~^ represents the cut 
we introduce to ensure there is no redshift dependent selection in the sam- 
ple. Throughout this paper we only consider the radio properties of BCGs 
with an iRadio— ccn greater than this (the radio-loud BCGs). The dotted 
line reoresents the 1 mJv threshold of the FIRST cataloeue. 
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2.1 OWLS Comparison simulations 

For comparison with theory we compare the resuhs in this pa- 
per with simulat ions from the Overwhelmingly Large Sim ulations 
project (OWLS. ISchave et alJl2010l : lMcCarthv et alj|2010h that are 
ideal for this study as they include the relevant physics, and output 
values for cluster mass, Lx, Tx, galaxy stellar mass and photome- 
try. The simulated Tx we choose for comparison with our observa- 
tions is emission weighted rather than spectroscopic-like, although 
for OWLS the difference between these is small with a median off- 
set of 0.1 keV. The volume-limited nature of the simulations means 
that they do not contain the massive clusters found in the XCS 
sample, but there is a significant crossover regime between the two 
(10^^ < Msoo < lO^M©). The two OWLS simulations we use 
are both high resolution, hydrody namical, cosmologica l simula- 
tions th at include radiative cooling ( Wiersma et al ]|200i), star for- 
mation iSchave & Dalla Ve cchia 200S ), feedb ack from supemovae 
driven winds foalla Vecchia & Schavd 1200 Sl) and full chemody- 
namics Jwiersma et al. 2009), but only one of the simulations in- 
cludes feedback from supermassive black holes (AGN feedback, 
iBooth & Schav3[2009i) . While the efficiency of the AGN feedback 
was set to reproduce the normalisation of the M — a relation, 
the model was not tuned to repro duce any ot her observables. As 
demonstrated by .McCarthy et al. (201(] ll201 ll) . in the AGN feed- 
back simulation the low entropy gas, which would otherwise con- 
dense to form stars, is removed or heated at high redshift (z > 1) by 
energy input from AGN in typical L* galaxies. This needs to hap- 
pen at high redshift as the gravitational binding energies are easier 
to overcome in the lower mass progenitors of the z — groups. 
Subsequent to the ejection of the low entropy gas, the injection of 
further energy by the AGN of the central galaxy stops the remain- 
ing gas from over-cooling in the centre of the group and thus stops 
the gas from causing significant star formation. 

Figure [2] illustrates the effect of this AGN feedback on the 
Lx — Tx rela tion within the OWLS simulations, as was also 
demonstrated in iMcCarthv et"aLH201 0Y Here, the lines on the plot 
connect the same halos from the non-AGN feedback run to those 
with AGN feedback included. We see that when comparing the 
AGN and non-AGN runs the effect of the AGN, to first order, is 
to lower the Lx for a given Tx- This is because the simulated 
AGN are heating and removing the low entropy gas, lowering the 
central density and thus Lx which depends on the integral of the 
square of the gas density. The influence of the AGN is more pro- 
nounced in the lowest Tx, and therefore lowest mass, groups as 
these systems have shallower gravitational potential wells and are 
therefore less able to hold onto their central gas. Similar predic- 
tions fo£_the_L_x_—7x relati on have been m ade by other groups 
(e.g. lPuchwein et alj!i2 008: Fa bian et alj2010l) . 

As in iMcCarthv et alj j2010l 1201 ll) , groups/clusters are se- 
lected on the basis of halo mass. Haloes are identified using a 
friends-of-friends (FOF) algorithm and we select only those haloes 
with M200 > IO^^Mq. Self-gravitating subhaloes are identified 
in eac h FOF system using the SUB FIND algorithm of Dolag et al. 
(l2009l). which is mo dified version of that originally developed by 
ISDringeletai] ( l200ll) . The BCG is defined as the stellar component 
associated with the most massive subhalo in a FOF system. Since 
our observations do not include ICL, we only quote the total stellar 
mass/luminosity of the BCG within a radius < 30 kpc. We find that 
the AGN feedback in the OWLS simulation reduces the stellar mass 
of the BCG by a factor of ~ 20 to bring them into agreement with 
values derived from observation. As a consistency check, we com- 
pare our derived mass-to-light ratio (1.91) used in ^to obtain the 
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Figure 2. A plot showing the effect of AGN feedback on the X-ray lumi- 
nosity - temperature relation in the OWLS simulation (as demonstrated in 
IMcCarthv et alj2010X To help illustrate this point the temperatures here are 
cool-core corrected. Blue squares are the non-AGN feedback run and red 
stars are with AGN feedback included. The connecting lines show where 
each individual cluster is moved to by the feedback process. The low en- 
tropy gas, which would otherwise condense to form stars, is removed or 
heated at high redshift (2 > 1) by energy input from AGN in typical L* 
galaxies. The ejection of low entropy gas and the injection of energy by 
the AGN in the central galaxy stops the remaining gas from over-cooling in 
the centre of the group and thus stops the gas from causing significant star 
formation in the central galaxy at late times. 



stellar masses for our observed BCGs and find the median value of 
that in the simulations to be in excellent agreement at 1.95 ± 0.03. 
The OWLS output data, presented as comparison throughout, is at 

z = 0. 



3 RESULTS 

3.1 BCG X-ray and mass scaling relations 

The optical properties of BCGs are k nown to cor r elate with the X 
ray properties of their host clusters ('Edge"l99lVCollins & Mam 
1998; Lin & Mohr 2004; Popesso et al. 2007; Whilev et al. 200^ 
Stott et al. 2008; Brough et al.ll2008l ; iMittal et al. 2009). These re- 
lations have been studied for large samples of massive clusters but, 
owing to the wide mass range of the XCS, we can now probe them 
for a consistently analysed sample of groups and clusters, spanning 
over three decades in Lx- A summary of the scaling relations is 
presented in Table[T] 

Figure [T] demonstrates that there is no evidence for a strong 
correlation in either Lx or Li-SDSS with redshift due to either 
an X-ray or optical flux limit (which could perhaps be present if 
cluster confirmation relied on bright BCGs). Satisfied that any cor- 
relations will not be due to a selection effect, in Figure |3| (upper 
and middle panels) we plot the Li_sDSS of the BCG against the 
Lx and Tx of the host clusters. This is the first time such rela- 
tions have been produced for a large sample of BCGs over such a 
wide range in Lx and Tx- There is clearly a positive relationship 
between BCG ii-sDss and both Lx and Tx- We q uantify this 
by performing a linear regression, BCES bisector fit dlsobe et al.l 
ll99d ; lAkritas & Bersh adv 1996) to the logged data, accounting for 
the errors in both quantities and the intrinsic scatter, finding that: 
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log {Li- 



(0.44±0.04) log {E(zy^ Lx) + (11.36±0.03)(2) 



log (i>_sDss) = (0.99 ± 0.05) log {Tx) + (10.74 ± 0.03), (3) 

where Li_sDSS is in units of Lq, Lx is in units of lO^^ergs"^, 
Tx is in keV and i;(2) = [Q.^{1 ^ zf + Q.aY''^ . 

We compare these findings with the OWLS AGN feedback 
run, plotting its output in the panels of Figure [3] (crosses). In gen- 
eral the OWLS simulation covers a lower cluster mass {Lx and 
Tx) range than the observations due to the volume-limited nature 
of the simulations. As discussed in i]2. 11 the simulated BCG pho- 
tometry and stellar mass give an average mass-to-light ratio consis- 
tent with that of our observed photometry and SSP derived masses. 
The 30 kpc radius aperture used to extract the simulated magnitudes 
is appropriate for comparison with the SDSS mo del magnitudes a s 
BCGs typically have effective radii of 30-40 kpc jStott et al.l2"oTTh . 
The plots show that the OWLS simulations are in good agreement 
with both of the observed X-ray - BCG luminosity data, as where 
there is an overlap, the Li_sDSS values agree and if taken together 
the two samples appear to form a continuous distribution. To quan- 
tify this we also perform fits to the OWLS data which yield BCG 

T ^ r0.40±0.02 . rynn T ^ rp0.81±0.04 „ -^v, 

.t^i-OWLS OC i^x-OWLS ™" BCU ivi-OWLS « Jx-OWLS' ^™ 

the former in excellent agreement with its observational counter- 
part. 

The relations seen in Figure [3] {upper and middle panels) 
demonstrate that BCG stellar mass correlates with the total mass 
of the cluster. To quantify this, we study the derived BCG stellar 
masses and cluster i\/5oo values calculated in ij2]and compared to 
the OWLS output. In Figure [3] (/ower panel) we plot BCG stellar 
mass against the cluster i\f5oo value and perform fits to the obser- 
vations and simulations. For the observations we find that the BCG 
stellar mass: 

log (M.bcg) = (0.78 ± 0.06) log (M500) + (11.19 ± 0.06), (4) 

with M500 in units of lO^^M©. This is important as in a cos- 
mological context it gives a relation between the stellar mass of 
the central galaxy and its dark matter halo mass over the range 
10" < Msoo < IO'^Mq. 

For the OWLS simulations M,bcg-owls oc M^^^^^^q^'^o,, 
again in excellent agreement with observation. The offset be- 
tween the fit zero-points is most likely due to the uncertainties in- 
volved in both deriving physical quantities from observables such 
as Li-SDSS and Tx and the difficulty in simulating realistic baryon 
physics. 

We compare our B CG cluster ma s s scali ng relation to that 
found in other studies. iHaarsma et alj jZOlCj) find that studies 
with photometry calculated in fixed metric apertures typically ob- 
serve shallower relations between BCG luminosity/mass and clus- 
ter mass than those that use isophotal magnitudes, with those using 
Sersic profile derived m asses, as is the case in this study, found 
to be the steepest CMittal et alj|2009l ). However, we also note that 
the large fixed 30 kpc radius aperture from the OWLS simulation 
does not appear to give a shallower slope. In the comparisons that 
follow we assume a linear relationship between BCG luminosity 
and BCG mass for the studies that only provide the dependency of 
BCG luminosity on cluster mass. The exponent a = 0.78 ± 0.06 
in the relation M«bcg oc M^i^^^^j. found in this study is indeed 
similar to the a = 0.62 iz 0.05 found by extrapolating Sersic pro- 
files in lMittaletal] ( l2009i) (who also use a BCES fitting routine). 
For the is ophotal and fixed m etric aperture studies a is significantly 
less ( e.g. iLin & Mohill2004l find a = .26 ± 0.04. |Popesso et al.l 
20071 find a = 0.2 5, IWhilev et alj|2008l find a = 0.12 ± 0.03 and 
Brough et alj2008l find q = 0.11 ± 0.10). The reason for this may 



be that the Sersic-like photometry gives brighter total magnitudes 
and perhaps therefore raise the value of a. 

BCGs can be described as galaxies at the centres of their dark 
matter halos and thus we compare to mass scaling relations de- 
rive d from clus t ering and halo occupation modelling. For exam- 
ple, lYang et al.l ( l2005h find L. oc M{^^i^ where a = 0.25 for 
halo masses above lO^^/i~^M0 when studying the halo occupation 
stati stics of cent ral galaxies in groups and a = 0.28 inferred by 
Vale & Ostriked 12004) from matching the galaxy luminosity func- 
tion to a theoretical halo mass function. It is clear that our a is sig- 
nificantly higher than these alternatively derived values, however, 
Zheng et al. (2009) find a ~ 0.5 from the clustering of central 
luminous red galaxies (LRGs) at z = 0.3 which is closer to that 
found here and a more appropriate comparison sample. 



3.1.1 BCG colours 

We also investigate BCG colours to look for those that are espe- 
cially blue or red by seeing how the g — r colour evolves with red- 
shift. Blue colours may indicate star formation, perhaps linked to 
cool-core activity jO'Dea et alj|200^ . Instead, we find that BCGs 
are remarkably uniform in colour with an r.m.s scatter of 0.09 mag 
across the entire redshift range, which is consistent with the photo- 
metric uncertainty. The colour evolution is in excellent agreement 
with the Bruzual & Chariot (2003) SSP used throughout the paper 
(Figure |4]l demonstrating that the stellar populations in BCGs are 
consistent with being formed at 2 = 3. We note that there is no 
correlation between colour and any other BCG or cluster property 
discussed throughout this study. 



3.2 The effect of BCG mass on the Lx — Tx relation 

As discussed in the introduction, the Lx — Tx relationship is 
used to understand the physics of the ICM with, to first order, 
Lx correlating with the integral of the gas density squared and 
Tx correlating with the virial mass of the system. In Figure |5] 
we plot the Lx ~ Tx relation for the groups and clusters in 
the XCS-SDSS sample. Indicated on the plot are the most lu- 
minous (red stars) and the least luminous (blue squares) 30% of 
BCGs (37 clusters in each). A linear regression fit to the logged 
Lx — Tx data, acco unting for the errors in both quantities (see 
iMehrtens et al. ll20lTlfor error values) and the intrinsic scatter, is 
performed for these two subsamples and the whole sample assum- 
ing the form log {E{z)~^ Lx) = a log (Tx) + 6, with Lx in units 
of 10"** erg and Tx in keV. As in i|3. li the fitting routine used is 
the BCES bisector method with errors deriv ed from bootstrapping 
dlsobe et al.|[l990l : lAkritas & Bershadvlll996 1. We note that the re- 
sults are not affected if we instead use the BCES orthogonal method 
which gives fit parameters consistent to within the la error of the 
preferred bisector method. The XCS selection function is not taken 
into account in these fits; however, there is little redshift depen- 
dence in the Lx limit over the redshift range considered (Figure 
[T] upper panel) and so we assume selection effects will not have 
a significant impact on our results (although see i]4.2| for a discus- 
sion on Malmquist bias). For the whole XCS-SDSS sample these 
best-fit parameters are a = 2.67±0.19 and b — — 1.54±0.07 (see 
also Hilton et al. in prepQ. This is typical of group and cluster sam- 
ples with a > 2 indicating a break from self-similarity dMushotzkvl 



^ Hilton et al. (in prep) will present the XCS Lx 
entire sample and its evolution with redshift. 



Tx relation for the 
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Table 1. BCG - Cluster scaling relation summary of the form log Y = a log X + b. 



Y 


X 


a 


6 


BCG Li_SDSs(-f'©) 


Tx(keV) 


0.99±0.05 


10.74±0.03 


BCG Li_SDSs(-^^0) 


Lx(10**ergs-i) 


0.44±0.04 


11.36±0.03 


OWLS BCG L,(Lq) 


Tx(keV) 


0.81 ±0.04 


10.74±0.03 


OWLS BCG L^\Lq) 


Lx(104^ergs-i) 


0.40±0.02 


11.42±0.05 


M.bCgCMq) 


M500(1014Mq) 


0.78±0.06 


11.19±0.06 


OWLS A^.bcg{M0) 


MsooClQi^M©) 


0.76±0.04 


11.47±0.04 



1984'; 'Edge & Stewart" 1991'; 'MarkevitchI 1 19981 ; lAmaud & Evrar j 
1 999; . Pratt et al. 2009; Mittal et al. 2011,) . 

For the luminous BCG I/i-SDSS sub-sample these parameters 
are aiu^-ao = 3.69 ± 0.49 and 6ium-30 = -2.02 ± 0.24 and for 
the faint BCG Lj-SDSS sub-sample they are afaint-30 = 1-92 ± 
0.33 and bfaint-ao = -1-45 ± 0.11. The Lx - Tx relation for 
the most massive BCGs is steeper than for the lowest mass BCGs 
to a significance of 3a, with a cross-over at Tx ~ 2 keV. There is 
therefore good evidence that above Tx ~ 2 keV, the most massive 
BCGs live in clusters with a higher Lx value for a given Tx than 
those that host the least massive BCGs, with the fits to the two 
populations crossing at Tx ~ 2keV. Clusters that host the least 
massive BCGs seem to follow a relation in agreement with the self- 
similar, a — 2, case. These results are summarised in Table|2] 

This steepening of th e Lx — Tx relat ion is also found in 
other sub-populations, with iMittal et al.l ^20091) finding strong cool- 
core c lusters have Lx oc y^-^sio is and lMagliocchetti & Briiggenl 
( l2007h finding Lx oc r* o=^°-^ for those containing BCGs with 
extended radio emission. Both these properties may be associated 
with clusters that contain the most massive BCGs, as we discuss in 

El 

For a theoretical point of view, we plot the Lx — Tx relation 
for the OWLS simulation and find a qualitatively similar situation, 
albeit over a different range in Tx/cluster mass owing to the lim- 
ited volume of the simulation (see Figure [SJ. With a and h repre- 
senting the same fit parameters as above, aowLS = 1-88 ± 0.16 
and &OWLS = -1-70 ± 0.07, aowLS-ium-so = 2.54 ± 0.21 
and bowLS-ium-30 = -1.45 ± 0.06, and aowLS-faint-ao = 
1.24 ± 0.13 and 6owLS-faint-30 = -1-98 ± 0.08. Although the 
slopes are systematically shallower, these fits still show the same 
relative behaviour, analogous to the observations. 

3.3 Dynamical state of the cluster 

One way of testing whether a cluster is relaxed, is to investigate the 
projected offset of the BCG from the X-ray centroid of the cluster 
in units of 7?5oo, the idea being that in the least disturbed, most viri- 
alised systems one would expect both the peak in the gas density 
and the BCG to be co- located at the centre of mass of the cluster 
jSanderson etal]|2009h . To measure this offset correctly we note 
that t he XCS and SDSS astrometry are matched ( Mehrt ens et al.l 
l201lh and that the optical centroiding of the BCG from SDSS is 
known to ~ 1 — 2 kpc. However, the XCS is a serendipitous sur- 
vey and as such many of the sources are off-axis and affected by the 
strongly varying point spread function (PSF) over the XMM New- 
ton field of view, which has to be corrected for fclovd -Davies et al.l 
l20Tlh . We therefore expect the main source of error will be the cen- 
troiding of the extended X-ray emission found by the XCS XAPA 
software. Because the XCS is serendipitous we can make use of the 
fact that a number of our sources are multiply imaged, some even 



being the on-axis main target in one exposure but off-axis in oth- 
ers. There are 13 such multiply imaged clusters in the redshift range 
of our sample and we report no correlation between off-axis angle 
and centroid offset, so this is clearly not an issue. However, there is 
a scatter in the centroid positions of the same clusters when mea- 
sured from different exposures, which has a median offset of 3.3 
arcsec, corresponding to an average offset for the sample of 1 1 kpc 
or~ 0.01i?50o. 

The median offset between the BCG and the X-ray centroid 
for the XCS-SDSS sample presented here is found to be 0.030 ± 
0.010 i?5oo- Figure [7] shows that there is no significant correlation 
between this offset and the luminosity of the BCG (I/i-SDSS oc 
BCGoffset"'' °^='=° ''^). This is surprising as one might expect that 
BCGs at the centre of clusters would have had more opportunity 
to accrete mass at the centre of the potential well, although, it is 
interesting to note that there are few low luminosity BCGs with 
small (< 0.01 -Rsoo) offsets from the X-ray centroid. 

For co mparison with other w ork, when we analyse the tabu- 
lated data of lHaarsma et d]j2010f), we find no correlation between 
BCG offset and BCG luminosity for their sample of X-ray lumi- 
nous clusters. They find that for 90% of their clusters the BCGs lie 
within 0.035 of the -R500 with all cool-core clusters within this dis- 
tance. However, for our sample this number is only 54% with 90% 
found within 0.10 -Rsoo- The reason for this difference ma y relate 
to the relative flux limits of the samples, with the Haarsma et al.l 
( I2OI0I) sampl e containing o nly m assive Mdustcr > lO^^M©, high 
Lx systems. iLin & Mohd ( |2004|) find that 95% of BCGs in their 
sample lie within 500 kpc of the X-ray centroid and we find this 
number to be 99% but this is to be expected as our clusters probe 
down to the group scale where 500 kpc is comparable to the virial 
radius. 

In Figure [8] (the equivalent of Figure [5j we plot the 
Lx - Tx relation for the XCS-SDSS sample. Fits of the form 
log (i5(2)~^ Lx) = alog(rx) + b, with Lx in units of 
lO^^ergs^^ and Tx in keV, are performed to the most offset 
30% of BCGs, with fldist-ao = 2.19 ± 0.19 and 6dist-30 = 
-1.30 ± 0.19 and least offset, with a„car-30 = 3.29 ± 0.51 and 

6„ear-30 = -1.72 ±0.18. 

In a situation mirroring that of the stellar mass of the BCG, the 
Lx — Tx relation for the least offset BCGs appears steeper than for 
the most offset BCGs to a significance of 2.1(7. There is therefore 
some evidence that above Tx ~ 2 keV the least offset BCGs live in 
clusters with a higher Lx value for a given Tx than those that host 
the most offset BCGs, again the fits to the two populations cross 
at Tx ~ 2keV. Clusters that host the most offset BCGs seem to 
follow a relation in agreement with the self-similar, a = 2, case. 
These results are summarised in Table[2l 

We speculate that the influence of the dynamical state is con- 
nected with strong cool-core clusters as they a re also found to 
have BCGs co-located with their X-ray centroids jSanderson et al.l 
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Figure 3. For the BCG - cluster scaling relation plots presented here, the 
observed and simulated data are represented as filled points and crosses re- 
spectively, and fits to the former and latter are solid and dashed lines. The 
radio-loud BCGs are plotted as open diamonds. The error bars for the X-ray 
data can be found in Figure lAll We note that the middle and lower panels 
are not independent - see ^for a description of how M500 was calculated. 
Upper: BCG SDSS i-band luminosity plotted against cluster/group bolo- 
metric X-ray luminosity for the XCS-SDSS sample. Middle: BCG SDSS 
i-band luminosity plotted against cluster/group X-ray temperature for the 
XCS-SDSS sample. Lower: BCG stellar mass plotted against cluster mass. 
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Figure 4. BCG SDSS g - r col our plotted against redshift for the XCS- 
SDSS sample (filled points). A iBruzual & CharloJ j20O3l) simple stellar 
population (SSP) with a Chabrier IMF, a formation redshift, z j = 3 and so- 
lar metallicity is plotted as the dashed line. The radio-loud BCGs are plotted 
as open diamonds. 

I2OO9I : iHaarsma et al]|2010l: IZhang et al.ll201lh and have a steeper 
Lx — Tx relation than the non-cool-core systems, with Lx oc 
y3.33±o.i5 j-Qj. cool-cores and Lx oc t^«±0 2i for 

non-cool-cores, in agreement with our findings for BCG offset 
jMittal et a i]i2or ij). The XCS is a serendipitous survey, it is based 
on cluster images taken over a wide range of off-axis angles (i.e. 
offsets from the instrument aim point) and so has degraded spatial 
resolution compared to a cluster survey based on targeted XMM 
follow-up (the PSF is strongly off-axis dependent). Therefore, we 
are unable to determine which of the 123 clusters in our study have 
cool cores, using spatial fitting. However, it seems a good assump- 
tion that a significant fraction of the l east offset BCGs will b e in 
cool-core clusters, giv en the findings of lSanderson et al.l ( l2009h and 
iHaarsma et"ai] ( l2010() . 

A further signature of the dynamical state of the cluster is the 
luminosity gap between the BCG and t he next brightest gal axy in 
the cluster (e.g. ISmith et al.ll201(]l see iHarrison et al.ll2012l for a 
detailed study of the luminosity gap in the XCS clusters). This lu- 
minosity dominance of the BCG over its neighbours should provide 
clues as to the merger history of the cluster, as a cluster with more 
than one 'BCG' may have undergone a recent merger, whereas one 
that is dominant over its satellites may be in a more relaxed system. 
We perform an analysis of the most dominant and least dominant 
30% of BCGs in the Harrison et al. sample within the context of 
the Lx ~ Tx relation, as above. There is a hint of an effect in that 
the most dominant BCGs have a steeper relation but it is not sta- 
tistically significant. With the fit parameters a and h as above, we 
find for the most dominant 30% of BCGs adom-30 = 3.72 ± 1.30 
and 6dom-30 = —1-90 ± 0.44 and for the least dominant this is 

Clnon-dom-30 = 2.71 ± 0.47 and 6non-dom-30 = -1.64 ± 0.15. 

3.4 Radio-loud BCGs 

As described in ^ we study the radio properties of the BCGs 
in our sample, considering 'radio-loud' BCGs to be those with 
^Radio-ccn > 2 X lO^^WHz"^, due to the flux limited na- 
ture of the FIRST survey (see Figure [T] lower panel). The frac- 
tion of radio-loud BCGs in our sample is 0.24 ± 0.05, which 
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Figure 5. The X-ray luminosity plotted against X-ray temperature. Red stars are the top 30% of clusters in terms of BCG SDSS i-band luminosity and blue 
squares are the bottom 30%. The dashed line is fit to the clusters with the most luminous BCGs, the dotted line is fit to those with the least luminous and the 
solid line is fit to all (see also Hilton et al. in prep). The luminous population is found to have a steeper slope than the faint population. The error bars for the 
X-ray data can be found in Figure IaT] 
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Figure 6. The X-ray luminosity plotted against X-ray temperature for the 
OWLS simulation. Red stars are the top 30% of clusters in terms of BCG 
i-band luminosity and blue squares are the bottom 30%. The dashed line is 
fit to clusters with the most luminous BCGs, the dotted line is fit to those 
with the least luminous BCGs and the solid line is fit to all. 



Figure 7. BCG SDSS i-band luminosity plotted against BCG offset from 
X-ray centroid. A fit to these data is represented with a dashed line which 
shows no significant congelation. The radio-loud BCGs are plotted as open 
diamonds. 



is similar to the fr action found in the X-ray selected sample of 
Lin & Mohj (l2007h (0.33) and in the optically selected sample of 
Best et al ](l2007|) (0.21±0.02 for clusters with velocity dispersions 
cr„ > 500 kms~^). As noted in <j2]these two similar samples define 



radio-loud at a fainter limit of iRadio 2i 1 x 10^^ W Hz~ but in the 
high stellar mass regime this fractio n is unaffected b y such a small 
change in the radio luminosity limit teest et alj2005l) . therefore we 
do not expect this to significantly affect any comparisons. In Figure 
|9]we plot the distribution of radio luminosities of the sample which, 



10 J. p. Stott et al. 




Figure 8. The X-ray luminosity plotted against X-ray temperature. Red stars are the clusters containing the 30% of BCGs least spatially offset from their 
hosts' X-ray centroids and blue squares are the 30% most offset. The dashed line is fit to those with the least offset BCGs, the dotted line is fit to those with 
the most offset BCGs and the solid line is fit to all. The co-located population is found to have a steeper slope than the most offset population. The error bars 
for the X-ray data can be found in Figure IaTI 



due to the flux limit, peaks around I/Radio-tot = 10 WHz 
with a maximum value of LRadio-tot = 1-2 x lO^^WHz"^. 

There is no evidence for a correlation between the total radio 
luminosity of the BCG and its optical luminosity or the mass of the 
cluster (Figure [Tot. One might expect to find a correlation between 
BCG I/i_sDSS and radio luminosity as stellar mass is thought to 
correlate with black hole mass via the Af — a relation and the accre- 
tion rate could thus be higher. In other studies of AGN in BCGs the 
maximum radio luminosity for a given stellar mass is a function of 
mass, but the relat ion for the entire sample has a large scatter, due 
to the duty-cycle jLin & Mohr 2007). In Figure [TTI (upper panel) 
we plot the fraction of radio-loud BCGs against BCG stellar mass 
but find only a weak dependence. We note that lBest et alj ( |2007|) 
find the BCG radio-loud fraction /rl oc ° when lower BCG 
masses are included and then a flattening off at high mass. So we 
can say that for our sample, BCG radio-loud fraction does not cor- 
relate wi th stellar mass bu t this is still consistent with the measure- 
ments of iBest et al.l ( |2007|) . because of the high stellar mass range 
considered. If we compare the radio-loud fractions above and be- 
low the median stellar mass of the sample (2.3 x lO^^Mo) then the 
most massive BCGs have /rl = 0.27 ± 0.07 and the least massive 
BCGs have a radio-loud fraction of 0.21 ± 0.06 so again there is 
no significant difference, demonstrating that stellar mass has very 
little influence on radio-loud fraction at these high stellar masses. 

There is also no correlation between the total radio luminosity 
of the BCG and its distance from the X-ray centroid. This is perhaps 
surprising as one might expect BCGs in relaxed clusters, being lo- 
cated closest to the centre of the cooling gas reservoir, would have 



more opportunity for fuel to be effectively fed to the black hole, but 
the duty-cycle may again be responsible for this result. However, 
in Fi gure [TTI (m iddle panel) we plot radio-loud fraction against dis- 
tance from the cluster X-ray centroid and find an excess towards the 
centre, suggesting a higher fraction of AGN hosting BCGs, in line 
with the hypothesis above. The fraction of radio-loud BCGs with an 
offset less than the median offset of 0.037?500 ./rl = 0.37 ± 0.09 
compared to 0.14±0.05 above this, a difference with a significance 
of 2.3a which indicates there is some connection between the two. 

Figure [12] shows the Lx — Tx relationship for the XCS- 
SDSS-FIRST sample and we find a fit to this relation of the form 
\og{E{z)-^ Lx) = a log(Tx) + h, where a = 2.73 ± 0.20 
and b = -1.53 ± 0.08 in agreement with the fit to the XCS- 
SDSS sample in Figure [5] For the radio-loud sample we find that 
QRL = 2.91 ± 0.45 and bRL = — 1.50 ± 0.21 in agreement with 
that of the entire sample. These results are summarised in Table|2] 
We include a dividing line at Tx — 2.1 keV which corresponds to 
the median Tx of the sample. 

The fraction of radio-loud BCGs in clusters with Tx < 
2.1 keV is /rl = 0.10 ± 0.04 and for Tx > 2.1 keV it is /rl = 
0.38 ± 0.09, a difference of 3a, demonstrating that cluster mass is 
an indicator of BCG radio-loudness (see Figure [TTj lo wer panel). 
We note that from scaling relations such as that in IXue & Wul 
( 2000) this value of Tx = 2. 1 keV (from the Tx -M500 scaling re- 
lation used throughout this paper, Afsoo = 1.4 X 1O^*M0) roughly 
corresponds to ay — 500kms~^ , the value of cluster velocity dis- 
persion used bv lBest etal.l ( l2007l) to separate their sample into high 
and low mass clusters and the M200 = 1-6 x 10^'* division used 
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Figure 9. The radio luminosity distribution for tlie radio-loud sample 

(iRadio-ccn > 2 X 1023WHz-l). 



bv lLm&Mohil(l2007h . The latter also find a significant increase in 
BCG radio-loud fraction with cluster mass as /rl = 0.36±0.03 in 
high mass clusters and /rl ~ 0.13 ± 0.05 in lo w mass systems, in 
excellent agreement with our findings. However. iBest et al] j2007h 
do not find a significant difference between the radio-loud fractions 
in their high- and low-mass cluster samples. The most probable rea- 
son for this discrepancy is that the Best et al. (2007) sample is op- 
tically rather than X-ray selected. This may mean that their sample 
includes a number of clusters which are perhaps under-luminous in 
X-rays for their optical richness and as such do not make it into our 
sample. Alternatively, as their clusters are selected on optical rich- 
ness, in order to measure a velocity dispersion, they are potentially 
only selecting the richest of the low mass systems which are per- 
haps more likely to host AGN. We note that we find no evolution 
in the AGN fraction with redshift, as for the low redshift half of 
our sample z < 0.2 /rl = 0.18 ± 0.06 whereas at z > 0.2 this 
number becomes /rl = 0.31 ± 0.08, a difference of only 1.3cr. 
Table[3]summarises the radio-loud fraction dependencies. 

To crudely transform these radio-lou d fractions in t o estim ates 
of AGN life-times using the method of iLin & Mohj i2007h . we 
can say that we know there is compelling evidence that BCGs 
are fully formed hy z — 1.5 iStott et al] I2OIQ) making them 
6.9 Gyr old at the median redshift of the survey, z — 0.2. This 
may be conservative as the BCG stellar populations formed at 
z > 3 (|Stot t et al. 20102, although potentially in separate sub-units 
jPe Lucia & Blaizotll2007h and the AGN fraction in clusters, and 
perha ps therefore duty-cycle, increases with redshift l lMartini et al.l 
l2009h . The low mass, Tx < 2keV (M500 < lO^M©, < 
500 kms~^) clusters therefore have life-time estimates of 0.7 Gyr 
and the high mass Tx > 2keV (M500 > 10^'' M©, cr„ > 
500kms~^) clusters have life-times of 2.7 Gyr. And for the entire 
population the life-time is 1.7 Gyr. 



4 DISCUSSION 

4.1 BCG AGN energy input into ICM 

To study the influence of the AGN on the ICM, we can compare 
rough estimates of the energy output of the AGN with the ther- 
mal energy of the ICM, following a similar prescription to that de- 
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Figure 10. BCG radio luminosity plotted against BCG SDSS i-band lumi- 
nosity (upper) and cluster mass (lower), both showing a lack of significant 
correlation. 



scribed in 'Birza n et al.1 ( |2004) and lLin&Moh3 j2007h . We see no 
significant correlation between BCG radio luminosity and cluster 
mass (Figure [Tot so, with the caveat of the lower duty-cycle (see 
^3.41 . we may expect that AGN outbursts will be more influen- 
tial on group scales. The mechanical power of the radio jet which 
causes 'bubbles' in the ICM is given by 400!^r7l/Radio-tot, where 
u is the frequency (1.4GHz), rj is the efficiency coefficient and 400 
is a factor that comes from the assumptio n that the gas conta ined in 
the outflows is relativistic and adiabatic (' Birzan et al.ll2004) . This 
power can be multiplied by the life-time of the AGN to give the 
integrated energy input over its history, .Eagn. For the average 
radio-loud AGN in our sample with median I/Radio-tot = 2.1 x 
10^'*WHz-i and a life-time of 1.7 Gyr, Eagn = 6.377 x lO'^^J. 
The thermal energy of the ICM is given by: 



-EiCM = (3/2)A;BTx/iCMM5oo//imp, 



(5) 



where /icm is the ICM mass fraction, the mean molecular weight 
/i — 0.59 and rrip is the proton mass. For the median cluster in our 
sample M500 = 1.4 x lO^^M© which corresponds to an /icm ^ 
0.09 of lGonzalez et all ( l2007h giving -BicM = 8.6 x 10 J. 

The ratio of the energy input of the AGN to the thermal energy 
in the whole system for an average cluster in our sample is therefore 
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Table 2. Lx — Tx relation summary of the form log Lx = o, log Tx + b, with Lx in units of 10**erg s~^, and Tx in keV. 



Sample 


a 


b 


XCS-SDSS ALL 


2.67±0.19 


-1.54±0.07 


XCS-SDSSTjf < 2.1 keV 


2.47±0.52 


-1.39±0.08 


XCS-SDSS 30% most luminous BCGs 


3.69±0.49 


-2.02±0.24 


XCS-SDSS 30% least luminous BCGs 


L92±0.33 


-1.45±0.11 


XCS-SDSS 30% most offset BCGs 


2.19±0.19 


-1.30±0.19 


XCS-SDSS 30% least offset BCGs 


3.29±0.51 


-1.72±0.18 


XCS-HARRISON 30% most dominant BCGs 


3.72±1.30 


-1.90±0.44 


XCS-HARRISON 30% least dominant BCGs 


2.71±0.47 


-1.64±0.15 


XCS-SDSS-FIRST ALL 


2.73±0.20 


-1.53±0.08 


XCS-SDSS-FIRST radio-loud 


2.91±0.45 


-1.50±0.21 


OWLS ALL 


1.88±0.16 


-1.70±0.07 


OWLS 30% most luminous BCGs 


2.54±0.21 


-1.45±0.06 


OWLS 30% least luminous BCGs 


1.24±0.13 


-1.98±0.08 



Table 3. BCG radio-loud fractions above and below median(X) 



X median(X) /rl < median(X) /rl > median(X) 

BCG stellar mass (Mq) 2.3 x lO" 0.21±0.06 0.27±0.07 

BCG offset (fisoo) 0.03 0.37±0.09 0.14±0.05 

Cluster Tx(keV) 2.1 0.10±0.04 O.38±0.09 



^?AGN/i?iCM ~ 0.007r; which means that for a ~ 1 X IO^Mq 
cluster the AGN has little effect on the extended ICM. However, 
this value can be a factor of ^ 2.5 higher in the centres of clusters 
(10% of virial radius) and the val ue of rj can also vary from 1/30 
- 20 and may change with time jBirzan et al. l l2004l : lLin&Mohil 
l2007h . The energy input of the AGN can potentially be a significant 
fraction of the thermal energy of the ICM at the centre of such a 
cluster. 

For a typical massive, M500 = 1 x IO^^Mq, cluster contain- 
ing a BCG with the median Laadio-tot = 2.1 x lO^^WKz"^ of 
our sample (as there is no evidence for an I/Radio-tot — M500 de- 
pendency. Figure llOt but with the enhanced life-time of our high 
cluster mass sample, found in ii3.4l to be 2.7 Gyr, and with a corre- 
sponding Tx ~ 8.0 keV and /icm = 0.13 ^Gonzalez et"d]|2007h 
we have -Bagn /Eicm ~ O.OOOSry. This is a factor of ~ 20 smaller 
than the AGN effect on the median mass cluster stated above. 

Alternatively, for a typical low mass, A/500 = 1 x 10^'^ M©, 
group again containing a BCG with the median I/Radio-tot of 
our sample but with the lower life-time of our low cluster mass 
sample found in ^3.4\ to be 0.7 Gyr and with a corre sponding 
Tx ~ 0.3 keV and /icm = 0.07 dConzalez et"ai]|2007h we have 
Sagn/Sicm = 0.4r;. This is comparable to the thermal energy of 
the entire ICM and thus the AGN can have a major influence on the 
gas physics of the group. 

The results of this over-simplified analysis are plotted for il- 
lustration in Figure [T3] From this plot we can see that a central 
radio source with the same radio luminosity in all clusters has sig- 
nificantly more effect on group scales than in massive clusters, with 
the AGN able to inject a significant fraction of the thermal energy of 
the ICM in the low est mass systems. A similar trend is also found 
by iMa et"al] ( I2OI ih when studying the energetics of AGN over a 
range of redshifts. 

If we consider the transition from AGN feedback to cooling 
dominated clusters in the Lx — Tx relations, discussed in ^3.2\ 



and i]3.3l that appears to take place at Tx ~ 2keV (see ^4.21 
in the context of Figure [13] then an average value of 77 — 5.5 
takes -Eagn/-Eicm to > 0.04 at the corresponding mass ~ 1 — 
2 X lO^'^ Mfl. If the value is boosted by a factor of ~ 2.5 in clus- 
ter cores jLin & Mohill2007D then Sagn/Bicm > 0.1 at masses 
below this and AGN feedback can start to contribute significantly 
to the core energetics, with Sagn/Z^icm > 1-0 in the cores of 
Afsoo < 4 X lO^^M© (Tx < 0.8 keV) systems which could re- 
move the central gas entirely. This simple energetic argument may 
explain the transition from cooling to AGN feedback dominance 
seen in i]3.2| and i]3.3| and discussed in i]4.2| 



4.2 Time averaged effect of AGN 

The scaling relations between the BCG stellar properties and the 
X-ray properties of the ICM demonstrate that central galaxy mass 
correlates with the mass of its host dark matter halo over the range 
10^^ < M500 < IO^'^Mq The steep relation we find is similar to 
that found bv lMittal et ail l l2009l) who use an almost identical def- 
inition of magnitude and the same fitting method. The agreement 
between our observed BCG - cluster scaling rel ations and those 
from the OWLS simulation with AGN feedback jMcCarthv et al.l 
I20T0I1 may mean that AGN feedback, invoked to regulate gas on 
both the galaxy and cluster scales, is indeed the dominant baryonic 
feedback process at work in dense environments. We note that the 
OWLS simulations indicate that ejection of low entropy gas at high 
redshift from L* galaxies is key to producing the dramatic change 
in the Lx — Tx relation seen in Figure [2l with the cooling of the 
gas at late times regulated by feedback from the BCG AGN. 

The lack of a significant correlation between BCG stellar mass 
and its distance from the X-ray cent roid adds weight to the obser- 
vation that B CGs form early 2; > 1 jColUns et alj2009l : lstott et al.l 
I2OIOI I20T l|) and that subsequent disturbances to the cluster, per- 
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Figure 11. Upper: The fraction of radio-loud BCGs per stellar mass bin for 
the XCS-SDSS-FIRST sample. The dashed line represents the congelation 
and subsequent flattening from lBest et al.l )2007fl . which is in good agree- 
ment with our findings. Middle: The fraction of radio-loud BCGs per bin 
in BCG offset from the X-ray centroid of the cluster for the XCS-SDSS- 
FIRST sample, which indicates that BCG radio-loud fraction increases 
when the BCG is co-located with the peak in the ICM surface brightness. 
Lower: The fraction of radio-loud BCGs per cluster mass bin for the XCS- 
SDSS-FIRST sample. The BCG radio-loud fraction increases significantly 
with cluster mass. 



haps via cluster-cluster merging, thus have little influence on BCG 
mass. 

The interrelation between the BCG properties and the ICM is 
highlighted in the fits to the Lx — Tx relation for different sub- 
samples. Fits to clusters which host the most and least massive 
BCGs are discrepant to a significance of 3a. The most massive 
BCGs live in clusters with a steepened relation, such that above 
Tx ~ 2keV (M500 ~ 1 X lO^'^M©) they reside in clusters with 
relatively high Lx compared to those with the lowest mass BCGs 
and below this they occupy the same region of the Lx — Tx plot. 
This is also seen in the OWLS simulations, albeit over a differ- 
ent cluster mass range, further demonstrating that simulations that 
model AGN feedback are broadly able to recreate observable prop- 
erties. The offset of the BCG from the X-ray centroid of the clus- 
ter has a similar effect on the Lx — Tx relation, with the BCGs 
co-located with their clusters' X-ray centroid seemingly having a 
steeper slope than their more offset counterparts. 

From the literature, other populations that have steepened 
Lx — Tx relation slopes are those with extended AGN radio emis- 
sion and those with a st rong cool-core iMagliocchetti & Briiggenl 
l2007l ; lMittaletal.ll201ll) . In a time averaged sense we expect the 
most massive BCGs in our sample to harbour the most power- 
ful AGN, even if they are not currently accreting. This is be- 
cause BCG mass correlates with black hole mass through the 
.Ai^BH — A/suigc relation and black hole mass governs the ac- 
cretion rate which is coupled to the AGN energy output (e.g. 
iBooth & Schavd 12010) . Having a BCG co-located with the X- 
ray centroid is also a ssociated with AGN activity and the pres- 
ence of a co oling core /Sanderso n et al.l2009l : lHaarsma et al.l2010l . 
iMittal et al . 201 1 ; Zhang et al. 201 1|), confirmed for this sample by 
the decrease in BCG radio-loud fraction with distance from X-ray 
peak (Figure[TT] middle panel). 

As Lx is a proxy for the density of the cooling gas, this steep- 
ening of the Lx~Tx relation means that in the cluster mass regime 
(Tx > 2 keV) the most massive BCGs, or those co-located with 
the cluster X-ray centroid, live in systems where gas cooling domi- 
nates, whereas in the group mass regime they live in systems where 
AGN feedback dominates. To explain this, we note that we showed 
in i]3.1| that BCG mass depends on cluster mass through a sub- 
unity exponent and that radio luminosity does not correlate with 
cluster mass ( i]3.4| l. Therefore, the effect of AGN feedback will not 
be self-similar In i]4.1| we see this in action as, assuming the same 
radio luininosity of AGN at each cluster mass, we find that its en- 
ergy becomes increasingly significant with decreasing cluster mass, 
despite the lower duty-cycle, simply because the gravitational po- 
tential well is shallower and the gas mass is lower in groups, mak- 
ing it easier to disturb the ICM and raise its temperature. BCGs 
co-located with the cluster X-ray centroid, and therefore with the 
peak in ICM density, will have a more effective supply of gas to 
feed their central black hole and as such the energetics argument is 
the same as that for the most massive BCGs. 

The clusters with the lowest mass BCGs, or those most offset 
from their X-ray centroid, have Lx — Tx relations that are con- 
sistent with the self-similar situation of Lx oc T\. So in clusters 
where there is little evidence that cool-cores and non-gravitational 
effects such as AGN feedback have an influence on the ICM, the 
X-ray observables show that this is indeed the case. 

The apparent lack of a departure from self-similarity in the low 
BCG mass/large BCG offset systems may give insight into their 
formation or recent merger history. Recent cluster-cluster mergers 
may have disrupted the potential, removing the opportunity for a 
cool-core to form, so the BCGs will look under-luminous as they 
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Figure 12. X-ray luminosity plotted against X-ray temperature for the XCS-SDSS-FIRST sample. Large black points are radio-loud (X/Radio— cen > 2 X 
lO^f W Hz^ ^ ) whereas small black dots are below this threshold. The dashed line is a fit to the radio-loud sample and the solid line to the whole sample. The 
veilical dotted line represent the median Tx of the sample. The error bars for the X-ray data can be found in Figure IaTI 



originate from one of the merging clusters and thus do not follow 
the average BCG mass-cluster mass trend from i ]3.1l A recent dis- 
turbance would also explain the offset of the BCG from the X-ray 
gas centroid. However, we find no correlation between BCG mass 
and BCG offset in i]3.3l Also, when analysing the effect of the lumi- 
nosity gap of the BCG down to the next brightest members, there 
is no evidence that clusters with several 'BCGs' follow different 
Lx — Tx relations. Alternatively, ignoring cluster-cluster merg- 
ing, perhaps these are systems where a strong mechanical or heat- 
ing mechanism removed the low entropy ga s at an early epoch (e.g . 
AGN feedback from L* galaxies at 2 > l. lMcCarthv etalj|201lh 
and this has been regulated ever since by AGN feedback from the 
BCG, making both the ICM and BCG under-luminous, in the mas- 
sive {Tx > 2 keV) clusters, due to the lack of cooling gas available 
for star formation. 

We note here that a fit to the group Lx — Tx relation {Tx < 
2.1 keV), gives a slope consistent with that of the whole population 
{Lx PC T^'^''^^'^^) in co ntrast with the steepened slope found by 
iHelsdon & PonmaiJ ( 12000 ) {Lx <x T^'^).We speculate that this 
difference is because that sample is optically selected, unlike the X- 
ray selected sample presented here. This lack of slope change on the 
group scale can be reconciled with the above results as, even though 
the gravitational potential is easier to overcome in groups, strong 
AGN feedback still only occurs in systems where it is required to 
stop gas over-cooling. 

A caveat is that our results are likely affected by selection bi- 
ases, although our sample has an Lx limit reasonably constant with 
2. If a lower portion of the group Lx range is missing at low Tx be- 
low this limit, then we may expect a shallower slope in this regime 



(i.e. see Figure 5 of lAllen et al.ll20l"ll) . which would particularly 
affect the Lx — Tx relation of those with low mass BCGs. We at- 
tempt to test this with the volume limited OWLS simulation as this 
is qualitatively similar to the observed data. By cutting the data at a 
Lx = 3 X lO^^ergs"^ to remove the 30% lowest Lx groups and 
refitting the Lx — Tx relation for the whole sample and the lumi- 
nous and faint BCG sub-samples. As one might expect the slopes 
all become shallower but are still consistent with the values from 
i] 2. II within Icr except for the slope of the faint BCG sample which 
changes by ^ 1.5(t and thus does not affect the conclusions. From 
this we conclude that the faint BCG sample will be affected most 
by selection bias, although it is difficult to estimate the magnitude 
of the effect, particularly with the complicated selection function 
due to the serendipitous nature of the XCS. 

4.3 Instantaneous effect of AGN 

The Lx — Tx relation fit to the radio-loud AGN is in agreement 
with that of the whole sample, although we do not have the statis- 
tics to go further and say whether this is associated with any of the 
optically selected sub-samples described in i]4.2| We find the result 
that AGN radio luminosity does not correlate with either the stel- 
lar mass of the BCG, its offset from the X-ray centroid, or cluster 
mass. As discussed in the results section, this is most likely due 
to the AGN duty-cycle which means that AGN are switched on and 
accreting intermittently and therefore any correlations of this nature 
are washed out. To understand how AGN relate to their host galaxy 
and c luster, the AGN fraction is a more useful measure teest et aU 
l2007h . 
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Figure 13. This plot is provided as an illustration of the energet ics derived 
from our simplified version of the model taken from lBirzan et al.l i2004D and 
iLin & MoiiJ )2007l) . The black points are the ratio of integrated AGN en- 
ergy to the thermal energy of the ICM plotted against cluster mass, assum- 
ing the same radio luminosity (iaadio— tot = 2.1 X lO'^^W Hz~^) across 
all masses, rj is the efficiency coefficient and the dashed line is a fit to the 
points. We speculate that when the energy injection from the AGN reaches 
10% of the thermal energy in the core, it can begin to have a significant in- 
fluence on the ICM and assume this value is reached at Tx ~ 2 keV which 
corresponds to the crossover between the steepened and self-similar-like 
relations in ^3.2l and i]3.3l We therefore estimate that r] ~ 5.5 on average 
so that, when in combination wit h a ~ 2.5 X heating enhancement in the 
central region iUn & Mohill2007l) . Eagn/Eicm > 0.1 at < 2 keV. 
Therefore the lower horizontal dotted line represents an estimate of where 
^AGN / ^ICM > 0. 1 in cluster cores. The upper horizontal dotted line rep- 
resents an estimate of where Sagn/^'ICM > 1.0 in cluster cores, which 
may lead to a significant disruption of the ICM at Tx < 1 keV. The shaded 
regions represent where ICM cooling dominates and where AGN feedback 
becomes influential, for this simple treatment. 



We find that the AGN fraction for the BCGs in our sample is 
only weakly dependent on BCG mass for the limited range of high 
galaxy masses we observe. One might expect to see a correlation 
as black hole accretio n rate is thou ght to be correlated with the 
mass of the black hole ( iBondill 19521) and thus the mass of the host 
galaxy through th e Mbh — Msuige relation iMagorrian et al.ll9"98h . 
iBest et al.1 j2007) do find a correlation with stellar mass but this 
is for a sample that encompasses BCGs down to a mass of 4 x 
10^" M0 whereas only 8% of our sample have stellar mass less than 
10^^ M0. The lack of correlation may indicate that in the cluster 
environment where all BCGs are massive, the mass of the black 
hole is not what is driving the AGN activity. 

The AGN fraction correlates most strongly with the cluster 
Tx and therefore cluster mass (and similarly with Lx)- The AGN 
fraction in clusters below the median halo mass (Tx ~ 2 keV, 
~ IO^'^Mq) being lower by 3a than those in clusters with a halo 
mass greater than this. An enhancement in radio-loud fraction is 
also found for BCGs within the projected median distance from the 
X-ray centroid compared to those outside this radius (0.03i?5oo)- 
Both of these phenomena relate to the amount of cooling gas avail- 
able as Lx correlates with gas density so as in H4.2\ this suggests 
that the black hole needs an environment with significant amount 
of cooling gas available to fuel it. The fact that the small number 
of radio-loud BCGs in the lowest mass (Tx < 1 keV) clusters are 
found in those with average to relatively high Lx may back this 



claim (see Figure[72l(. As radio-loudness is an indicator of ongoing 
AGN feedback we may expect these clusters to perhaps lower in 
Lx or increase in Tx over the duty-cycle of the feedback process. 

From the above relations between AGN fraction and both 
galaxy and cluster properties it seems that the key parameters that 
govern the presence of AGN in BCGs are primarily the cluster 
mass/Tx/Lx and to a lesser extent the BCG offset from the cluster 
X-ray centroid, but not BCG mass. A picture is therefore emerging 
that the supermassive black holes at the centres of BCGs in clus- 
ter cores know more about their host cluster than they do about 
their host g alaxy. Whil e this is consistent w ith some simulations 
teooth & S chave 201ol: ICaspari et al.ll201 ll) and results from the 
clustering of radio-loud, luminous red galaxies ( Wake et al. 200^, 
it is perhaps puzzling as it means that the sub pc scale is governed 
by the 100s kpc scale and not the kpc scale. We conclude that the 
connection must be due to the amount of fuel available and thus 
what drives AGN activity in the central galaxies of clusters is the 
availability of a large, cooling gas supply that can be effectively fed 
to the central black hole. This explains the cool-core AGN connec- 
tion and the break from self-similarity for these sub-populations. 
How this cooling gas on kpc scales is funnelled down into the ac- 
cretion disk of the black hole remains an open question. 

The correlations found in this paper imply that we can use 
complementary information about the component galaxies to char- 
acterise sub-samples in current and future cluster surveys such as 
the XCS which compromise on X-ray integration time in favour of 
large area coverage, in order to improve constraints on cosmologi- 
cal parameters. 



5 SUMMARY 

The XCS first data release jMehrtens et al.|[201lh is key to this 
study as it provides an ideal sample of groups and clusters that span 
a large range in both Lx and Tx, which is essential to probe the 
effects of AGN feedback on the ICM in a large, statistically sig- 
nificant, sample. The results lead us to conclude that Tx ~ 2 keV 
(Msoo ~ IO^^Mq, Lx ~ 10*^ergs-\ cr„ ~ SOOkms"^) is an 
important mass scale which delineates the boundary between the 
cooling dominated systems above this mass and the AGN feedback 
dominated systems below. We propose that this is an excellent em- 
pirical demarcation between what is referred to as a 'group' and 
what is a 'clu ster' in agreement with p revious, perhaps ad hoc, def- 
initions (e.g. iBower & Baloghl 120041) and similar to the break in 
the entrop y-at-a-radius-enclosin g-a-given-mass vs. temperature re- 
lations of iPonman et al. I j2003h . Recent results indicate that AGN 
cavity pow er exceeds the luminosity of the cooling I CM at a similar 
mass scale jO'SulIivan et al.ll201 iMcitti et alj201 ih . 
We summarise our key findings as follows: 

1. The Lx — Tx relation for the clusters with the most massive 
BCGs is significantly steeper than for those with the least massive 
BCGs, with a similar story for those least and most offset from the 
peak in the ICM emission. This must relate to both AGN and cool- 
core phenomena as similarly steep relations are found for these sub- 
samples in the literature. 

2. Clusters without a massive, centrally located BCG instead fol- 
low self-similarity, confirming that it is the complex interplay be- 
tween cooling gas and the central AGN that leads to the departure 
from simple gravitational heating. We speculate that the self-similar 
clusters may have a different feedback, formation or recent merger 
history which stops gas from cooling onto the BCG's black hole. 
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3. The steepened and self-similar Lx — Tx relations described 
in points 1. and 2. cross at Tx = 2 keV and thus we conclude that 
gas cooling dominates above this mass scale and AGN feedback 
dominates below. 

4. We derive BCG-cluster scaling relations over the mass range 
10" < Msoo < 10" M©. The A/.bcg — Afsoo relation has an 
exponent < 1, implying that BCGs do not grow at the same rate as 
their host dark matter halos. The X-ray scaling relations can be suc- 
cessfully reproduced by cosmological hydrodynamical simulations 
that include AGN feedback, such as OWLS. 

5. A quarter of the BCGs in our sample are radio-loud, in agree- 
ment with similar studies, but we find a strong dependence on en- 
vironment, in that those in clusters with Tx > 2 keV are 4 times 
more likely to host a radio-loud AGN than those with Tx < 2 keV. 
The radio-loud fraction also increases for BCGs closest to the peak 
in the ICM emission but the stellar mass of the BCG has little effect. 
This indicates that it is the availability of a fuel supply of cooling 
gas from the halo environment that governs AGN activity in the 
most massive galaxies. 

6. From a crude calculation of the energetics to explain points 
1., 2. and 3., it is easier for AGN to affect the group scale systems 
than the cluster scale, which may be as a direct result of point 4., 
the sub-unity exponent in the relationship between BCG mass (and 
therefore black hole mass) and the mass of its host cluster com- 
bined with the lack of correlation between BCG radio luminosity 
and cluster mass. 

7. The XCS sample shows no evidence for the group scale Lx — 
Tx slope being different from that of the cluster scale. This can be 
reconciled with the above results as, even though the gravitational 
potential is easier to overcome in groups, strong AGN feedback still 
only occurs in systems where it is required to stop gas over-cooling. 
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APPENDIX A: ADDITIONAL FIGURES 



BCG AGN feedback and the ICM 




Figure Al. The X-ray luminosity plotted against X-ray temperature for the XCS-SDSS sample to illustrate their associated errors. 



